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DIFFRACTION  EFFECTS  IN  DIRECTED  RADIATION  BEAMS 


I.  Introduction  and  Summary 

Diffraction  is  a  fundamental  characteristic,  of  all  wave  fields,  be  it 
photons,  electrons,  etc.  The  effect  of  diffraction  is  typically  manifested 
when  an  obstacle  is  placed  in  the  path  of  a  beam.  On  an  observation  screen 
some  distance  away  from  the  obstacle  one  observes  a  rather  complicated 
modulation  of  the  time-average  intensity  in  the  vicinity  of  the  boundary 
separating  the  illuminated  region  from  the  geometrical  shadow  cast  by  the 
obstacle. 

In  many  applications  it  would  be  highly  desirable  to  propagate  a  beam 
over  a  long  distance  without  an  appreciable  drop  in  the  intensity.  As  an 
example  we  cite  the  possibility  of  accelerating  particles  to  ultra-high 
energies  by  utilizing  high -power  laser  beams.  Although  the  accelerating 
gradient  in  many  of  these  schemes  is  extremely  large,  the  actual  distance 
over  which  the  particle  and  laser  beams  maintain  an  appreciable  overlap  is 
very  limited.  The  overlap  is  reduced  due  to  the  diffraction  of  the  laser 
beam  and  as  a  result  the  net  gain  in  the  particle  energy  is  limited. 

With  the  advent  of  high-power  lasers  and  microwave  sources, 
diffraction  of  radiation  beams  with  finite  transverse  dimensions  has  turned 
into  a  problem  cf  special  importance.  As  an  example,  consider  laser 
radiation  of  frequency  w  emanating  from  a  cavity  oscillating  in  the 
fundamental  transverse  Gaussian  mode.  How  far  will  this  beam  propagate  in 
a  turbulence- free  atmosphere?  More  to  the  point,  how  fast,  is  the  fall-off 
in  the  intensity  of  this  laser  beam? 

i 

The  answer  to  this  question  is  well-known.*  The  scale  length  for  the 
fall-off  in  intensity  is  given  by  t  fie  Rayleigh  range,  defined  by 
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where  w^  is  the  minimum  spot  size,  or  radius,  of  the  beam,  and  X  =  2nc/«  is 
the  wavelength.  The  minimum  spot  size  wq  is  also  known  as  the  waist  of  the 
radiation  beam.  The  fall-off  of  the  beam  intensity  as  it  propagates  in 
space  is  a  consequence  of  the  fact  that  initially  the  beam  was  constrained 
to  a  finite  waist,  wq.  Diffraction  then  causes  the  beam  to  spread  in  the 
lateral  direction  and,  from  energy  conservation,  the  intensity  must  drop 
off  correspondingly.  In  the  limit  of  an  infinitely  wide  beam,  w^  ->  <*>,  the 
Rayleigh  range  is  infinite,  there  is  no  diffraction  and  the  intensity  is 
constant . 

A  natural  way  to  propagate  a  beam  over  long  distances  is  to  increase 
the  Rayleigh  range  by  employing  a  wider  beam  or  shorter  wavelength 
radiation.  Clearly  the  width  of  the  beam  is  limited  by  the  energy  source 
available  for  pumping  the  lasing  medium,  and  short  wavelength  lasers 
(x-rays  and  beyond)  are  not  presently  available.  As  a  result,  over  the 
past  several  years  there  has  been  an  upsurge  in  research  on  such 
fundamental  topics  as  propagation  and  diffraction  properties  of  radiation 
beams.  (See  Ref.  2  for  an  earlier  discussion.)  Briefly,  the  question 
being  asked  is:  "Can  diffraction  be  overcome?"  The  following  is  a  summary 
of  our  review  of  diffractionless  and  other  directed  radiation  beams. 

i)  Electromagnetic  Missiles  (Section  IV) 

Experiments  indicate  the  possibility  of  generating  wave  packets  with  a 
broad  frequency  spectrum.  The  high-frequency  end  of  the  spectrum 
determines  the  furthest  distance  the  missile  can  propagate,  in  complete 
accord  with  our  understanding  of  diffraction. 

ii)  Bessel  Beams  (Section  V) 

A  Bessel  beam  is  a  particulai,  monochromatic  solution  of  the  wave 
equation.  Bessel  beams  propagate  no  fintber  than  Caussian  beams  oi  plane 
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waves  with  the  same  transverse  dimensions  and,  contrary  to  previous 
assertions,  Bessel  beams  are  not  "resistant  to  the  diffractive  spreading 
commonly  associated  with  all  wave  propagation". 

iii)  Electromagnetic  Directed  Energy  Pulse  Trains  (Section  VI) 

These  are  particular,  broad-band  solutions  of  the  wave  equation.  We 
show  that  the  experiment  and  the  numerical  studies  of  these  pulses  are 
consistent  with  conventional  diffraction  theory  and,  contrary  to  previous 
assertions,  these  pulses  do  not  "defeat  diffraction". 

iv)  Electromagnetic  Bullets  (Section  VII) 

Electromagnetic  bullets  are  solutions  of  the  wave  equation  which  are 
confined  to  a  finite  region  of  space  in  the  wave-zone.  The  ultimate  goal 
of  the  research  has  been  to  determine  the  source  function  which  leads  to  a 
prescribed  form  for  the  bullet  in  the  wave-zone.  Although  the  mathematical 
framework  for  this  has  been  established,  no  concrete  example  has  appeared 
in  the  literature. 

Sections  II  and  III  begin  with  a  review  of  basic  diffraction  theory 
and  our  findings  and  conclusions  are  summarized  in  Section  VIII. 
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1 1 .  Elec  tromagnetic  Wave  Pi  f  f  rat:  t  iort 

Consider  the  radiation  beam  from  a  cavity  of  radios  d.  The  wave 
vector  is  given  by  k ( ( e z ,  corresponding  to  propagation  predominantly  in  the 
z  direction,  and  the  magnitude  of  the  spread  in  the  wave  vector  in  the 
transverse  direction  is  denoted  by  Ak^ ,  with  >>  Ak^.  The  angular  spread 
of  the  radiation  relative  to  the  z  axis  is 

0  =  Ak^/k(| .  (2) 

On  an  observation  screen  at  a  distance  z,  the  radius  of  the  illuminated 
region  is  given  by 

w  =  d  +  9z .  ( 3 ) 

The  first  term  on  the  right-hand  side  of  this  expression  indicates  the 
width  of  the  region  illuminated  according  to  geometrical  optics.  Beyond 
this  lies  the  region  of  the  geometrical  shadow,  and  the  second  term  in 
Eq.  (3)  indicates  the  extent  to  which  this  region  is  illuminated  due  to 
diffraction  of  light.  The  distance  2  over  which  the  angular  spread  leads 
to  a  fall-off  in  the  intensity  is  given  by  d  +  9Z  =  2d,  or 

Z  =  d/6.  (4a) 

The  distance  Z  may  be  regarded  as  the  scale-length  for  diffractive 
spreading  of  the  beam. 

As  a  first  example,  suppose  the  transverse  distribution  of  intensity 
in  the  beam  is  uniform.  This  is  the  case  when  plane  waves  are  apei lured. 

If  the  radius  of  the  aperture  i:  ,  tom  a  undamental  result  of  Fourier 
analysis,  Ak^d  -  1.  The  angulat  spread  is  therefore  given  by 
0  =  A/2nd,  where  A  =  2  n  /  k  ? ,  if;  ihe  wavelength.  For  this  intensity 
distribution  one  thus  finds 
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(4b) 


Zp  =  2nd2/X. 

For  the  case  when  the  transverse  intensity  distribution  is  a  Gaussian, 

2  2 

exp(-r  '/vo>,  of  width  w  ,  ve  have  Ak^  =  1/vq,  an^  angular  spread  of  the 

beam  is  on  the  order  of  9  =  X/2itw  .  In  this  case  d  =■•  w  and  hence 

o  o 

ZG  =  2rtw2/\,  (4c) 

which  is  twice  the  Rayleigh  range  Zp  defined  in  Eq .  (1). 

Clearly,  diffraction  is  simply  the  physical  manifestation  of  the  well- 
known  result  of  Fourier  analysis  relating  the  spreads  in  wave  vector  space 
with  the  corresponding  widths  in  real  space,  Ak.jAx..  =  1  for  i  =  1,2,3.  As 
a  result,  Eq.  (4a)  expresses  a  fundamental  relation  which  we  shall  make  use 
of  repeatedly  in  order  to  interpret  the  results  of  theory  and  experiment  on 
so-called  diffractionless  radiation  beams. 
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Ill .  Diffraction  Zones  (Huygens'  Principle) 


According  to  Huygens'  principle  each  point  on  a  given  wavefront  acts 

as  a  source  of  secondary  wavelets.  The  field  at  a  point  P  is  given  by  the 

sum  over  the  wavelets.  If  u(r)e  lM)*  is  the  amplitude  on  an  aperture,  an 

3  A 

approximate  solution  of  the  scalar  wave  equation  at  P  is  given  by  ’ 

i  /  „  \  c-w-i  -iwt  P  n-l  iooR/c  /c. 

^p(r,t)  =  (iX)  e  j  dS'  u(r')  Re,  (5) 

aperture 

2  2  21/2 

where  R  =  [(x-x')  +  (y-y')  +  z  ]  is  the  distance  between  the  area 

element  dS'  on  the  aperture  and  the  point  P,  as  shown  in  Fig.  1. 

In  the  Fresnel  approximation  the  binomial  expansion  of  R  may  be  used 
to  simplify  Eq.  (5)  to 


*  (iXz) 


-1  iw(z-ct)/c 


Jd 


dS'  u ( x ' ,y ' )  e 
aperture 


.  (a)  f  x '  2  +  y ' 

1c  {  2z 


(6) 


For  plane  waves  incident  on  an  aper  *  ure  with  linear  dimension  d,  there 
are  two  physically  interesting  limits  for  approximating  Eq.  (6). 
i )  Fraunhofer  Diffraction  (Tar-Field  or  Wave- Zone  Region) 

If 


z  >>  d2/X,  (7) 

one  may  neglect  the  quadratic  terms  in  the  exponent  of  Eq .  (6)  and  the 
wavelets  from  the  entire  wavefront  at  the  aperture  contribute  to  the  field 
at  P.  In  the  Fraunhofer  region  ip  is  simply  the  Fourier  transform  of  the 
amp  1  i  t  udo  a  t  t  he  d  i  f  t  i  ac  t  i  ng  aper  tint  . 
ii)  Fresnel  Diffraction  (Near  Field  Regina) 

1  ii  t  he  other  1  i  m  i  t  , 

7.  <  d ‘  / X ,  (Hi 
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it  is  necessary  to  retain  the  quadratic  terms  in  the  exponent  of  Eq .  (6) 

and  the  wavelets  from  a  limited  portion  of  the  wavefront  at  the  aperture 

make  the  dominant  contribution  to  the  field  at  P.  In  this  case,  the 

integration  in  Eq .  (6)  may  be  taken  to  be  over  the  entire  z  =  0  plane. 

For  plane  waves  incident  on  a  circular  aperture  of  radius  d,  Fig.  2, 

making  use  of  Eq.  (5),  the  exact  field  on  the  axis  of  symmetry  is  given  by 

2  2  1/2 

4/p  exp(-iwt) {exp(  i«z/c)  -  exp[iw(d  +  z  )  '/c]},  and  the  intensity 

★ 

1  “  ^P^'p  is 

I  11  1  -  cos  {«((d2  +  z2)1/2  -  z]/c}.  (9) 

Figure  3  is  a  schematic  plot  of  the  intensity  function,  Eq .  (9),  indicating 

in  particular  the  transition  between  the  Fresnel  and  the  Fraunhofer 

regions.  Note  that  the  intensity  drops  off  precipitously  beyond 
2 

z  -  2nd  /X,  consistent  with  the  scale-length  defined  by  Ea .  (4b). 

Ve  proceed  now  to  examine  the  research  on  new  solutions  of  the  wave 

b  -20 

equation,  with  particular  emphasis  on  their  diffraction  properties. 


IV .  Electromagnetic,  Missiles 


i )  Theory 

Consider  first  the  case  of  a  field,  termed  a  issile",  which  falls 

off  more  slowly  than  the  usual  1/R  law.  The  inventive  step  is  the  use  of  a 

broad  frequency  spectrum.  Depending  on  the  spectrum,  the  fall-off  with  R 

7  8 

may  be  as  slow  as  desired.  ’ 

To  appreciate  the  nature  of  this  field,  note  that  for  an  arbitrary 
source  distribution  within  a  region  A  as  shown  in  Fig.  4,  the  energy 
delivered  to  a  screen  S,  integrated  over  all  time,  is 


e(S,R) 


|  dt  JdS  n  •  (ExB)/4 n, 

_<*>  screen 


where  n  is  a  unit  vector  normal  to  the  screen,  and  E  and  B  are  the  electric 
and  the  magnetic  field,  respectively.  For  a  source  with  a  bounded 

frequency  spectrum,  a  screen  of  fixed  area  S,  and  for  sufficiently  large  R, 

2  21 
f:(S,R)  ~  1/R  ,  according  to  well-known  results. 

The  current  density  for  the  eli  omagnetic  missile  described  in 

Ref.  7,  J(r,t)  -  &(z)f(t)ex>  r  <  d,  is  confined  to  a  disk  of  radius  d, 

?  2  1/2 

where  r  -  (x  •»  y  )  is  the  radial  coordinate  and  f(t)  is  a  given 
function  of  time.  It  A(oo)  and  J( w)  denote  the  Fourier  transforms  of  the 
vector  potential  and  the  current  density  respectively,  then 

A(to)  --  Jd  i  J(  w)oxp(  i  wR/c ) /c:R  is  a  solution  of  the  wave  equation.*"^  In  the 
present  case,  J(w)  -  S(z)f(w)e  ,  and  the  Vectoi  potential  on  the  axis  of 
symmetty  is  given  by 


A  ( u>) 


f'(w) 


i  (  t 


d  i  i  (  t  ‘ 


A 

V 


K 


Making  use  of  this  expression  for  A(«),  the  Poyntirig  flux  along  the  z  axis 
integrated  over  all  time,  U(z)  s  Jdt  e^ .  (ExB) /4n,  is  given  by 


-2  2 

C  Z[1  +2 ( Z  ^ 


2  -!/'> 
d/ )  1  " 


>1 


du>  f  ( w) 


cos 


rw -  ,  2 
{-1(2  + 


,2,1/2 

d  )  -z 


»)• 


(10) 


Note  the  resemblance  between  Eq.  (9)  and  the  integrand  of  Eq .  (10). 
Equation  (9)  is  for  a  monochromatic  field  and  is  based  on  Huygens-Fresnel 
principle  while  Eq.  (10)  is  obtained  from  a  rigorous  solution  of  the  full 
wave  equation. 

In  the  limit  z  -*  ®  in  Eq.  (10),  for  a  fixed  frequency  w, 

2  2  1/2  2 

cos  { (to/c )  ( ( z  +d  )  -z]}  ■+  cos(wd  /  2cz)  ->  1,  and  the  integrand  tends  to 

zero.  Referring  to  Fig.  3,  this  means  that  the  contribution  of  this 
frequency  lies  in  the  far- field  region  and  is  thus  negligible.  At  a  given 
large  z  we  can,  therefore,  write 

tr>  2  co  2 


12  J 


dw 


f(«) 


[1  -  cos(oxl2/2cz)]  =  ~ 


du> 


£(w) 


2  c  z  /  d ' 


2  c  z  /  d 


We  see  that  the  most  important  contribution  is  I i om  the  high  frequency  end 
of  the  spectrum,  foi  which  the  given  point  z  lies  in  the  near-field, 
Fresnel,  zone.  The  contributions  from  all  the  lower  frequency  components 
will  have  decayed  to  negligible  values  before  reaching  the  given  z.  While 
the  eventual  fall  off  of  any  frequency  component  is.  as  1/z,  the  fall  off  of 
the  time  integrated  Povnting  flux  for  tin  wave  packet  depend.’,  uii  how 
rapidly  |f  (w)  |  decay,  for  <•>  >  2  c  z  /  d  ‘  .  As  an  example,  corn:  idet  a  source 
with  a  ftequenry  tespon.se 

f  (  to)  “  [  1  *  (  ui-'  w  )  ’  |  (1  ^  * , 
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where  e  >  0.  For  this  spectrum  the  time-integrated  Poynting  flux  falls  off 
2  e 

as  U  “  1  / 7.  '  .  The  fall-off  can  thus  be  as  slow  as  desired  by  taking  the 

limit  e  -*  0.  The  limit,  however,  corresponds  to  a  frequency  spectrum  from 

a  source  with  infinite  energy.  The  function  f(t)  is  symmetric  with  respect 

to  t  ->  -t.  Evaluating  the  inverse  Fourier  transform  of  F(w),  one  finds 

that  f(t)  ~  t  ^  C~3 )  ^exp( -00  t)  for  t  >>  1/w  .  For  e  >  1/2  and  t  <<  1  /  to  , 

K  o  o  o 

f(t)  -»  constant.  However,  for  the  more  interesting  case  of  e  <  1/2,  f(t)  -> 

e-1/2 

t  when  t  <<  l/“0»  indicating  a  mild  singularity, 

i i )  Experiment  (Electromagnetic  Missiles) 

The  difficulties  involved  in  an  experimental  study  of  electromagnetic 
missiles  stem  from  the  need  to  generate  pulses  with  extremely  short  rise- 
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times  and  suitably-shaped  wavefronts.  An  antenna  was  used  to  generate  a 

"pure"  spherical  wave  which  formed  the  primary  pulse  and  the  field 

reflected  from  a  parabolic  dish  of  radius  2  ft  formed  the  secondary  pulse. 

The  pulses  were  detected  by  a  specially-designed  sensor.  The  primary  pulse 

2 

was  found  to  fall  off  as  1/z  ,  the  energy  decaying  by  1/16  when  the  sensor 
was  moved  from  4  ft  to  16  ft  from  the  source.  This  was  because  the 
antenna,  being  a  point  source,  generated  spherical  wavefronts.  The  pulse 
reflected  from  the  parabolic  dish  was  found  to  resemble  that  of  a  circular 
disk,  similar  to  that  studied  earlier  in  this  section.  Over  the  same 
distance,  the  energy  in  this  electromagnetic  missile  was  found  to  decay  by 
just  under  1/2.  Without  a  precise  knowledge  of  the  frequency  spectrum  it 
is  not  possible  to  make  a  quantitative  analysis  of  this  experiment.  Rough 
estimates  indicate  that  the  scale  length  ior  the  fall  off  of  the  intensity 
of  the  missile  is  indeed  compatible  with  the  diffraction  scale  length 
Z  ...  2nd' /A,  where  d  =2  ft  is  the  radius  of  the  reflecting  dish  and  a  is 
the  wavelength  for  the  highest  frequency  (10  GHz'  in  the  pulse. 
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These  preliminary  experimental  results  indicate  that  a  suitably 
tailored  pulse-shape  can  be  designed  to  have  an  energy-decay  rate 
essentially  limited  by  the  highest  frequencies  present  in  the  pulse 
generator,  in  complete  accordance  with  the  elementary  notions  of 
diffraction  of  light.  Propagation  of  a  composite  pulse  in  free  space  is  a 
dispersive  process.  As  the  beam  propagates,  the  lowest  frequency 
components  diffract  away  first. 


V. 


Bessel  Beams 


i )  Theory 

An  example  of  a  so-called  diffractionless  electromagnetic  beam  is  a 
Bessel  beam.  Ve  note  that  a  particular  solution  of  the  scalar  wave 
equation 


(V2  -  <f2  £_)  *  (r,t)  =  0, 


is 


vp  =  e 


2rt+e 

i(k|jZ  -  wt)  p  0  ik^(xcos9  +  ysin9) 


j  d6  A(9)e 


(ID 


2  2  2  2 

for  arbitrary  0q  and  A(0),  provided  <r'  =  c  (k(|  +  k^).  Here,  k|(  and  k^ 
denote  the  magnitudes  of  the  components  of  the  wave  vector  parallel  and 
orthogonal  to  the  z  axis,  respectively,  and  X  =  2.Jl/(kn  +  k.)  is  the 
wavelength.  Since  the  z-dependence  in  Eq .  (11)  is  separated  from  the  x- 
and  y-dependence,  the  solution  is  clearly  diffractionless  in  the  sense  that 
the  time-average  intensity  is  independent  of  z.  I"  fact,  the  intensity  is 
constant  for  all  7  and  all  t. 

Durniri  considers  the  case  where  A(Q)  =  1  (Ref.  10).  In  this  case, 
making  use  of  the  expa >sion  exp(itsinO)  =  E  J  ( C)exp(inQ) ,  where  J  is  the 
ordinary  Bessel  function  of  the  first  kind  of  order  n  (Ref.  22;,  Eq.  (11) 
simplifies  to  ^  -  2nJ  (k^r)expl  i(k|(z  -  art)  ) ,  where  1  ( y/  ,  y‘  )w  is  the 

radio!  variable. 


Making  use  of  the  properties  of  the  Bessel  function  (Ref.  22),  one  can 

2 

show  that  the  "energy"  content,  Jdr  1  J  '  (k  (  r ) ,  integrated  over  any 
transverse  period,  or  lobe,  is  approximately  the  same  as  that  in  the 


central  lobe.  This  point  will  be  important  in  our  interpretation  of  the 
diffractive  properties  of  Bessel  beams, 
i i )  Experiment  (Bessel  Beams ) 

A  Bessel  beam  has  an  infinite  number  of  lobes  and,  therefore,  has 
infinite  energy.  In  the  laboratory  an  approximation  to  this  ideal  beam  is 
realized  by  clipping  the  beam  beyond  a  certain  radius.  The  question  is, 
given  the  finite  transverse  size,  how  well  is  the  diffractionless  property 
preserved . 

To  answer  this  question  Durnin  et  a'l.  compared  the  propagation  of  a 
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clipped  Bessel  beam  with  a  Gaussian  beam.  The  full  width  at  half-maximum 

(FWHM)  of  the  Gaussian  was  taken  to  be  equal  to  the  FWHM  of  the  central 
lobe  of  the  Bessel  beam.  In  the  experiment  the  on-axis  intensity  of  each 
beam  was  measured  along  the  axis  of  symmetry.  The  Bessel  beam  was  claimed 
to  be  "resistant  to  the  diffractive  spreading  commonly  associated  with  all 
wave  propagation"  since  its  intensity  was  observed  to  remain  approximately 
constant  for  a  much  longer  distance  than  the  Gaussian  beam.  The  idea  of  a 
diffraction-free  beam  was  further  reinforced  by  using  a  geometrical  optics 
argument  to  obtain  a  formula  for  the  propagation  distance  of  central  lobe 
of  the  Bessel  beam. 

We  shall  now  reconsider  this  comparison.  The  wavelength  of  the 

o  2  2 

radiation  was  A  =  6328  A.  For  >he  Gaussian  beam.  exp(-r  /v£),  v  was  equal 

to  0.042  mm,  corresponding  to  a  FUHM  of  0.07  mm.  For  the  Bessel  beam, 

Jo(kjr),  was  equal  to  41  mm  \  corresponding  to  a  FWHM  for  the  central 

lobe  of  0.07  mm.  The  beams  were  aper lured  to  a  radius  d  -  3.5  mm.  The 

following  order -of -magnitude  discussion  is  based  on  Eqs.  (2)  (4);  a  more 

rigorous  analysis  is  presented  in  the  Appendix.  The  angular  spread  due  to 

the  natural  width  of  the  Gaussian  beam  is  6  A/2nw  and  Eq .  (3)  takes 

o 


i  he 


form  w  =  v  +  (X/2itv  )z,  where  the  first  term  in  this  expression  is  w  , 
o  o  o 

rather  than  d ,  since  the  energy  of  the  Gaussian  is  concentrated  in  the 

central  peak.  The  scale  length  for  diffraction  is  the  same  as  that  given 

y 

by  Eq.  (4c),,  namely  ^  2/rv^/X  =  1.75  cm.  The  natural  angular  spread  of 

the  Vessel  beam  is  6  =  kjX72n,  and  Eq.  (3)  takes  the  furm  w  =  d  + 

(k.^\/2n:)z,  where  the  first  term  represents  the  radius  of  the  aperture  since 
the  energy  in  each  lobe  is  approximately  the  same  and  they  all  affect  the 
propagation  of  the  Bessel  beam.  The  scale  length  for  diffraction  is, 
therefore,  given  by  d  +  (k^X/2n)Zg  =  2d,  or  Zg  =  2nd/k^X  «  85  cm,  which  is 
consistent  with  the  experimental  observation. 

In  the  transverse  plane  the  lobes  of  the  Bessel  beam  diffract  away 

sequentially  starting  with  the  outermost  one.  The  outermost  lobe  diffracts 

2  2 

in  a  distance  on  the  order  of  2ji  /Xlt^,  which  is  approximately  equal  to  Z^,. 
The  next  lobe  diffracts  away  after  a  distance  on  the  order  of  ZZ^.  This 
process  continues  until  the  central  lobe,  which  diffracts  away  after  a 
distance  ~  NZ^,  where  N  denotes  the  number  of  lobes  within  the  aperture. 

In  the  experiment  N  =  50,  implying  a  propagation  distance  on  the  order  of 
50Z„  for  the  central  lobe  of  the  Bessel  beam,  which  is  consistent  with  the 
measured  value.  Measurements  of  the  on -axis  intensity  obviously  fail  to 
reveal  the  gradual  deterioration  of  the  transverse  beam  profile,  but  the 
numerical  plots  in  Fig.  2  of  Ref.  10  are  consistent  with  this  scenario. 
Therefore,  the  Bessel  beam  is  not  "resistant  to  the  diffractive  spreading 
commonly  associated  with  all  wave  propagation".  Our  interpretation  points 
out  the  significance  of  each  successive  lobe  having  about  the  same  energy. 
The  centra]  lobe  persists  as  Jong  as  there  are  of  f -axis  1 obes  compensating 
for  its  energy  loss  and  hence  the  comparison  with  the  narrower  Gaussian 
beam  in  Ref .  1 t  is  of  little  significance. 
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Ve  note  that  utilizing  the  full  width  of  the  aperture  a  Gaussian  beam 
propagates  a  distance  on  the  order  of  NZ^;  i.e.,  N  times  further  than  the 
Bessel  beam.  Additionally,  by  appropriately  curving  the  wavefront,  nearly 
ail  the  power  of  the  Gaussian  beam  can  be  focussed  on  a  target  of  dimension 
w^  in  a  distance  Zg.  Hence,  for  this  purpose  a  Gaussian  beam  would  be 
significantly  better  than  the  Bessel  beam  employed  by  Durnin  et  al .  (See 
also  Ref.  12) 


IS 


VI.  Electromagnetic  Directed  Energy  Pulse  Trains 


i )  Theory 

Electromagnetic  directed  energy  pulse  trains  are  particular  solutions 

13 

of  Maxwell's  equations.  To  discuss  these,  we  make  the  change  of 
variables  £  =  z  -  ct  and  T  =  t,  and  transform  the  wave  equation 

(v2  -  c-2  32/ 3tz j  Y  =  0, 

into  the  form 


2  3 

c  3!,3t 


-2  2  21 
c  3  /3t  Y  =  0. 
y 


Making  the  assumption 

./i-  v  iwEyc 
V  =  K,  r ,  t)e  s  , 

leads  to  the  an  equation  for  the  complex  envelope  Yr 


+2io>e 


-2  3 


3x  +  c 


«r_ 

3t3~C 


-2 


32  1 
^  - 


0. 


(12) 


Here,  r  denotes  the  radial  variable  and  is  the  differential  operator  in 
the  plane  z  =  constant.  If  »K£tr,T)  varies  slowly  compared  to  the 
characteristic  scales  1/w  and  c/'w,  the  second  derivative  of  the  envelope 
function  may  be  neglected  anc  the  wave  equation  reduces  to 


(„2  .  -23),  n 

V,  +  2  lw  c  -  0. 


(13) 


Equation  (13)  is  an  extremely  useful  approximation  to  the  full  wave 
equation  in  a  vacuum.  Note  that  the  full  wave  operatoi  is  of  the 
hynerbolie  type,  where  ;  the  reduced  wave  operator  is  of  the  parabolic- 
type.  For  this  rt  ison,  Eq .  (13)  is  .sometimes  refer  red  to  as  the  parabolic 
approximation  to  tfie  wave  equation. 


lb 


A  particular  solution  of  Eq.  (13)  is  given  by 


-1  2  2 
v  -itan  ( t/ Tn )  -  ( 1 - i t/ td ) r  /w 

,  r,  °  K  K 

*  =  C  ~  e 


where  C  is  a  constant, 


1/2 


2.1 

v  *  wo  |^1  +  (t/tr)  J 


is  the  spot  size,  wq  is  the  waist  and 


2  2 
t_  =  ww  /2c  , 
R  o 


is  related  to  the  Rayleigh  range  Z 


R 


tow  /2c  by  t0  =  ZD/c. 
o  J  R  R 
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Ziolkowski  makes  use  of  the  variables  transformation 


£  =  z  -  c t ,  n  =  z  +  ct, 
in  the  wave  equation  to  reduce  it  to  the  form 

K  * 4  lk)T  *  °- 

Representing  Y  in  the  form 
Y  =  ^h,r)elwJ-/c, 

leads,  without  any  approximation,  to  an  equation  for  ip, 

K  *  41  “  Ifi  ]*  "  °' 

A  particular  solution  of  Eq .  (16)  is  given  by 


*  - 


C 


w  i  tan 

o 


(,VV 


( 1  i  q/  qR )  1 


(14a) 


(14b) 


(14c) 


(15) 


(16) 


(17a) 


where  C  is  a  constant, 


W  ---  W 


(3.7b) 


i  +  (n/^)7 


1/2 


and 


Hg  .  MW0/C, 


(17c) 


2 

is  related  to  the  Rayleigh  range  =  (tjvQ/2c  by 


R’ 


Some  remarks  on  the  solutions  in  Eqs.  (14)  and  (17)  are  in  order. 
First,  Eq.  (14)  is  a  solution  of  the  parabolic  approximation  to  the  full 
wave  equation.  On  the  other  hand,  Eq .  (17)  is  an  exact  solution  of  the 
full  equation.  Second,  there  is  a  factor-of-tvo  difference  between  the 
scale  length  cx^  in  Eq.  (14c)  and  the  scale  length  in  Eq.  (17c).  Third, 
the  solution  in  Eq.  (17)  has  infinite  energy.  Finally,  the  exact  solution 
in  Eq.  (15)  consists  of  a  pulse  traveling  to  the  left  which  is  modulated  by 
a  plane  wave  moving  to  the  right. 

To  examine  the  last  two  points,  Eqs.  (15)  and  (17)  may  be  combined  to 
form  a  fundamental  Gaussian  pulse  with  parameter  k  oo/c 


tk(rf z, t) 


ikq 

e 


-kr?7V 
e _ 

4niV"  ’ 


(18) 


where 


111 

V  A  ~  R  ’ 

A  =  z  +  £ } !  z  ,  R  =  f,  +  z 2  /  i , 

C)  o  o 

1 3  1 5 

ant,  z  is  a  constant.  To  conform  to  Ziolkovski '  s  example,  F.q .  (10) 

represents  a  pulse  traveling  to  the  tight  which  is  modulated  by  a  plane 
wave  moving  to  the  lefr.  With  an  appropr  iate  weight  lunction,  it  can  he 


IK 


shown  that  the  for  all  k  form  a  complete  set  of  basis  functions,  each 
with  total  energy  proportional  to  Jd'\|Yk|2  ->  «.  Just  as  in  the  case  of 
Fourier  synthesis  with  plane  waves,  a  general,  finite-energy  pulse  may  be 
obtain  ’  by  superposing  the  various  Y^  according  to  a  weight  function  F(k), 
that  is, 


00 


f  ( r ■ ,  z ,  t ) 


J  dk  Yk(r,z,t)  F(k) 
o 


■  4g(F-r  rtr  J dk  F<k)  e"ks-  <19a> 

o 

o 


where 


s  =  -if)  + 


(19b) 


Equation  (19a)  indicates  that  f(r,z,t)  is  proportional  to  the  Laplace 
transform  of  F(k). 


i  i )  Modified  Power-Spectrum  Pulse  ---  Numerical  Stu  dy 

Ziolkovski  has  examined  in  detail  the  pulse  corresponding  to  a 
modified  power-spectrum  (MPS): 


F(k)  -  ( 

Uni  g 


igjf  - 


r(«) 


e  -  a  ( (3k.  -  b) 


0  <  k  <  b/p 
k  >  b/p 


(20a) 

where  Ha)  is  the  Gamma  function  (Ref.  22)  and  a,  a,  and  0  are  arbitrary 


constants.  Upon  substituting  Eq -  (20a)  into  Eq .  (19)  one  obtains 


f  (  r  ,  ■£  ,  t  ) 


i  l 


In  /  g{  z  +  i  l.)  . .  . 

t.  £  ,  tbn/p  -j 

i  *  /  H(  7  t  i  £,)  1 1 1/  0 )  °  J 
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a  -+ 


(70b) 


The  real  part  of  this  function  defines  the  MPS  pulse.  The  radial  profile 
of  the  MPS  pulse  at  the  pulse  center,  £  =  0,  has  the  form 

,  ibp/3  -br2/gz 

f(r ,z)  =  - - ~ - - e  °.  (20c) 

Zo  (a  +  r2/gzo  -  .i0/3)a 

In  the  numerical  studies  the  pulse  vas  replicated  by  superposing  the 

fields  from  a  planar  array  of  discrete  points,  each  of  which  was  driven  by 

a  function  specified  by  the  MPS  form  on  some  z  =  constant  plane. ^  The 

10  -1  15 

parameters  were:  a  =  1.0  cm,  b  =  1.0x10  cm  ,  g  =  6.0x10  ",  •-= 

-  3 

1.667x10  ,  and  a  =  1.  The  spectrum  was  approximately  flat  up  to  290  GHz, 

becoming  negligible  beyond  15  TFz.  The  pulse  generated  in  this  manner  was 
propagated  forward  and  compared  with  the  exact  form  in  Eq.  (20b)  at  several 
locations  along  the  z  axis.  The  minimum  radius  of  the  array  required  to 
replicate  the  exact  pulse  form  at  1,  10,  100,  and  1,000  km  was  determined. 
From  Ziolkowski 's  results  we  estimate  the  corresponding  radii  of  the 
antenna  to  be  approximately  0.5,  5,  50,  and  500  m,  respectively. 

We  shall  examine  these  results  by  asking:  What  is  the  scale-length  for 

diffraction  of  the  MPS  pulse?  The  pulse  has  a  Gaussian  radial  profile,  as 

,  1/2 
indicated  in  Eq .  (20c),  with  a  width  w^  =  (gz  /b)  =  31.6  cm  and, 

therefore,  Ziolkowski  calculates  a  Rayleigh  range  m//A  =  0.21  km  for  the 

200  GHz  component . This,  however,  is  not  the  appropriate  scale- length 

for  diffraction  of  the  MPS  pulse.  The  correct  scale- length  is  given  by 

2nw  d/A,  where  the  antenna  dimension,  d,  always  exceeds  w  .  The  point  here 
o  o' 

is  that  the  Rayleigh  range  based  on  t he  waist  w  ,  as  calculated  by 

|  ^  -j  IR 

Ziolkowski,  ’  is  only  valid  at  the  pulse  center ,  C,  •  U.  Away  l  r  oin  the 
plane  {,  t)  the  effective  waist  increases,  as  indicated  by  F.q .  (20b),  and 
the  actual  difftaction  length  is,  theroloic,  longc-t  than  nv  ‘  A .  This 
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explains  why  the  pulse  propagates  further  than  the  Rayleigh  range  defined 

in  terms  of  v  .  To  calculate  the  actual  diffraction  length,  we  note  that 

the  perpendicular  wave  number  (k^)  spectrum  given  in  Ref.  15  indicates  that 

the  smallest  k,  is  on  the  order  of  1/w  .  Hence,  an  estimate  for  the 
1  o 

diffraction  angle  is  X/2nv^.  The  width  of  the  radiation  beam  given  by  Eq. 
(3)  can  be  written  as  w  =  d  +  fXZ2rtvo)z,  where  d  is  the  radius  of  the  array 
or  "antenna."  The  scale-length  for  diffraction  is  then  simply 


Zunf  =  2rtw  d/X. 
MPS  o 


(21) 


Note  the  similarity  between  the  diffraction  length  in  Eq.  (21)  and  the 
scale-length  for  diffraction  of  the  Bessel  beam,  Zg  =  2ttd/k^X,  derived  in 
Section  V,  subsection  ii).  The  resemblance  is  a  reflection  of  the  fact 
that  in  both  cases  the  pulse  energy  is  spread  over  the  entire  radius,  d,  of 
the  aperture,  which  is  much  larger  than  the  nominal  waist  of  the  beam,  wq. 

According  to  Eq.  (21)  the  larger  the  radius  of  the  array  is,  the 
longer  the  distance  of  propagation  of  the  pulse,  consistent  with  the 
numerical  results.  From  the  numerical  results,  the  ratio  ZMpg/d  is  equal 
to  2,000  which  is  the  same  as  that  given  by  Eq.  (2.1)  provided  the  frequency 
is  300  GHz.  Since  this  frequency  is  well  within  the  cutoff  of  the  pulse 

spectrum,  this  constitutes  a  persuasive  indication  that  the  MPS  pulse  does 
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not  'defeat  diffraction"  as  claimed  by  Ziolkowski. 

iii)  Modified  Power-Spectrum  Pulse  --  Experiment 

Ziolkowski  et  al.  have  performed  a  water-tank  experiment  to 

demonstrate  the  properties  of  a  MPS  acoustic:  pulse.  ^  The  pulse  was 

generated  by  a  6x6  chi’*  square  array.  Tire  MP.S  pulse  parameters  were 

1  4 

a  1 .0  m,  b  -  600.0  rn  ,  P  -  300.0,  z  4.5x10  in,  and  u  1  .  From 

o 

these  parameters  one  finds  that  the  pulse  width  v  is  equal  to  1.5  cm. 


The  experiment  indicated  that  a  Gaussian  pulse  with  an  initial  width  equal 

to  1.5  cm  suffered  a  greater  transverse  spreading  than  the  MPS  pulse. 

This  experiment  may  be  examined  in  the  light  of  the  discussion  leading 

to  Eq .  (21).  The  expression  in  Eq.  (21)  gives  the  scale-length  for  the 

fall-off  in  the  intensity  of  a  pulse  which  is  generated  by  an  array  (i.e., 

2 

antenna)  of  radius  d.  Since  the  square  array  is  6x6  cm  ,  we  take  the 

parameter  d  to  be  equal  to  3  cm.  Noting  that  the  speed  of  sound  in  water 
3 

is  1.5x10  m/s.  the  wavelength  of  the  dominant  frequency  in  the  pulse,  0.6 

MHz,  is  X  =  2.5  mm.  From  this,  the  actual  diffraction  scale-length  Zuri„  is 

1.1  m.  This  is  in  good  agreement  with  the  experimental  observation  that 

the  MPS  pulse  propagated  a  distance  of  1  m  without  significant  spreading. 

2 

Comparing  the  MPS  pulse  generated  by  a  6x6  cm  array  with  a  Gaussian  pulse 
having  a  waist  of  1.5  cm  is  inappropriate.  A  Gaussian  beam  with  spot  size 

equal  to  the  array  radius  used  in  the  experiment  would  prop. gate  a  distance 

2 

rtd  /X  =  Z^pc.;  i  •  £ •  as  far  as  the  MPS  pulse. 

Ve  note  from  Eq .  (21)  that,  in  general,  a  Gaussian  beam  with  an 

appropi lately  curved  wavefront  and  an  initial  spot  size  equal  to  the 

antenna  dimension  transfers  nearly  all  the  power  onto  a  target  of  dimension 

w  in  a  distance  on  the  order  of  /...nr.-  Such  a  Gaussian  beam,  therefore, 
o  MPS 

transfers  more  power  on  the  target  than  the  corresponding  MPS  pulse. 


VII.  Electromagnetic  Bullets 


In  this  section  we  shall  discuss  solutions  of  the  wave  equation  which 
are  confined  to  a  finite  region  of  space  in  tne  wave  zone  and  are  termed 
"electromagnetic  bullets".  Ve  consider  solutions  of  the  wave  equation 

(v2  -  c'2a2/at2)f(r,t)  =  -  p(r,t),  (22) 

where  the  source  term  p(r,t)  is  assumed  to  be  non-zero  for  a  finite  time 
interval  -T  <  t  <  T.  In  this  problem  there  are  two  cases  of  interest: 

Case  (a)  is  the  direct  source  problem  (initial  value  problem).  In 
this  case  the  solution  for  t  <  T  is  given  and  the  solution  for  t  >  T  is 
sought . 

Case  (b)  is  the  inverse  source  problem.  Here,  the  solution  of  the 
homogeneous  wave  equation  for  |t|  >  T  is  known  and  one  seeks  the  source 
term  appropriate  to  this  solution.  This  case  is  of  particular  interest 
since  it  would  enable  one  to  find  the  time-dependent  source  for  a 
prescribed  radiation  field. 

The  following  four  subsections  summarize  the  extensive  research  of 
Moses  and  Prosser  on  this  subject.  For  a  brief  description  of  the 
properties  of  a  bullet,  the  reader  is  referred  to  subsection  iv). 
i  )  Non  uniqueness  of  the  Inverse  Source  Problem 

In  this  subsection  we  indicate  the  reason  for  the  non  uniqueness  of 
the  inverse  source  problem. By  making  use  of  the  eigenf  ret  ions  of  the 
cut  i  operator  ,  t  he  electromagnetic  vector  tic  Id  wave  equation  may  tic*  solved 
along  essentially  the*  same*  lines  as  t  he  one  dimensional  pi  oh  1  cm.  Tire 
discussion  in  this  suhMn  tion  is,  therefore,  confined  to  the  one 
dimensional  wave*  equation  t  >  avoid  the  complications  of  multi  dimensional 


e I f  er t  s . 


Let  f  (x,t)  =  f(x,r)  denote  the  solution  of  Eq.  (22)  for  t  >  T,  ami 
£  (x,t)  =  f(x,t)  denote  the  solution  for  t  <  -I.  It  is  well-known  that  the 
solution  of  the  so  "-ce-free  initial  value  problem  for  t  >  T  in  terms  of  the 
values  of  the  function  f  (x,t)  at  r  =  T  and  the  "velocity"  (<)/3t)f  (x,t)  at 

•f  4- 

t  --  T  is  expressible  in  terms  of  a  propagator  G(x;t): 

f+(x,t)  -  Jdx'  [G(x-x';t-T)|^f+(x',t=T)  +  f+(x' , t=T)f^G(x-x' ; t-T) J . 

(23a) 

Similarly,  the  solution  of  the  source-free  final  value  problem  for  t  <  -T 
in  terms  of  f_(x,t)  at  t  =  -T  and  ( 8/ 9t ) f _(x , t )  at  t  =  -T  is  given  by 

f  (x,t)  =  dx'  JG(x-x' ;  t+T)|— £  (x',t.--T)  +  f  (x'  ,  t=-T  )f-G(x-x'  ;  t+T)  ] . 

-  J  0  t  —  -  O  L 

(23b) 

The  propagator  G  can  be  written  in  terms  of  the  Heaviside  function  q  as 
follows 

G(x;t)  =  |  sgn(t)  q(c"t2-  x2)  =  |  [q(x+ct)  -  q(x-ct)J.  (24) 

That  is,  G  may  be  expressed  as  the  difference  between  the  advanced  and  the 
retarded  Green  functions.  As  a  result,  f  (x,t. )  is  influenced  only  by 
points  x  at  t  ■-  T  which  lie  in  the  backward  light  cone  of  the  observation 
instant;  similarly,  f  (x,t)  is  influenced  only  by  points  x  at  t  -  -T  which 
lie  in  the  forward  light  cone  of  the  observation  instant,  as  indicated  in 
Fig.  5. 

Let  us  now  consider  the  effect  of  the  source  on  the  solution  of  the 
wave  equation.  \’e  define  two  auxiliary  functions, 

CO 

o(k, t )  (2n)  j  dx  p(x, t  )  e 


( 2  :> ) 


f  (k ,  cf,  t )  =  ?(k,txf-T)e  +  T) 


i  a  -iffclklt  f  ,  ...  ioc 

+  — irjkj  e  J  ^  p(k,  l  )  e 


fk.tn 


(26) 


where  a  -  +1  or  -1  distinguishes  the  two  directions  of  propagation  along 
the  x  axis.  Note  that  p(k, t )  is  simply  the  spatial  Fourier  transform  of 
p(x,t).  In  terms  of  the  two  auxiliary  functions,  it  is  simple  to  show  that 


cu 

f(x,t)  =  (2n)”i/2  Y,  J  dk  eikx, 

C  -«> 

is  a  solution  of  the  inhomogeneous  wave  equation  in  Eq .  (22).  The 
solutions  for  t  <  ~T  arid  for  t  >  T  are  then  given  by 


f+(x, t) 


.. 0  v~l/2  r< 

(2n)  u 


dk  e 


ikx  i cc  jk  j (±T  -  t ) 


f ( k, a, +T) .  (27) 


Equation  (27)  expresses  the  solution  of  the  wave  equation  in  terms  F 
evaluated  at  +T.  This  function  is  related  to  the  source  p  by  Eqs.  (25)  and 
(26).  Equation  (23),  on  the  other  hand,  expresses  the  same  solutions  in 
terms  of  f  and  (3/3t)f  evaluated  at  fT.  Hence,  one  would  expect  to  be  able 
to  relate  f(k,o,+T)  with  f  and  (9/3*:)f  evaluated  at  ±T.  Indeed,  the  the 
formula  connecting  F(k,or,T)  and  f  is 


C« 

F ( k ,  e , T )  =  (8n)  1/2  j  dx  [f4(x,i--T)  +  1  |f 7  flf  i  ( x ’  *  ~T )  ]  e  lkX>  (28) 

....  cl 


and  F(l;,o,  T)  is  obtained  in  terms  ol  f  by  Fourier  inversion  cf  F.q .  (27) 
evaluated  at  t  =  T.  Thus ,  in  t  he  direct  source  problem,  i  (x ,  i  )  in 
obtained  by  specifying  ei the i  f  ( x  ,  t  )  and  (3/3t)f(x,t)  evaluated  at  t  T 
oi  ' ( x , ~T )  and  p(x, t ) . 
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We  now  turn  to  the  inverse  source  problem.  Inverse  problems,  .in 
general,  have  been  the  subject  of  extensive  research  in  many  branches  of 
physics.  In  our  case  we  wish  to  determine  the  source  from  a  knowledge  of 
the  field  generated  by  that  source.  Supposing  that  f+(x,t)  and  T  are 
known,  we  can  determine  f(k,o,±T)  from  Eqs.  (27)  and  (28).  Letting  the 
upper  limit  of  integration  in  Eq,  (26)  equal  T,  it  appears  that  one  can 
then  obtain  the  temporal  Fourier  transform,  s(k.,w),  of  p(k,t),  where 

T 

s(k,«)  -  (2ji)-1/2  |  dt'  p(k,t')  eiwt'. 

•-T 


The  source  function  is  then  given  by 


p(x,  t) 


(2rc)  1 


co  cr 


dk  j  d«  s(k,<o)  eilcX 


Thus,  if  s ( k , w)  were  known  for  all  k  and  all  co,  wt  could  determine  p(x,t). 
However,  referring  to  Eq.  (26),  we  notice  that  s(k,co)  is  only  known  for 
co  *  ±  cjk],  which  is  not  sufficient  to  reconstruct  p(x,  t) .  This 
complication  is  intimately  connected  with  the  non- uniqueness  of  the  inverse 
problem. 

Moses  has  shown  that  specification  of.  the  time- dependence  of  the 

1 7 

source  is  sufficient  to  guarantee  a  unique  solution.  As  a  concrete 
example,  if  p(r,t)  =  p  (r)h  ( t )  -*  p  (r)h  ( t) ,  where  h  is  an  even  function 

of  t  and  h  is  an  odd  function  of  t,  and  both  essentially  arbitrary,  then  a 

o 

complete  solution  of  the  inverse  problem  for  p  (r)  and  p  (r)  is  possible. 


It  must  he 

s t  ressed , 

however,  that  this  as: 

umed  form  ?r 

>r  the  : 

:» o  *  i  ree 

function  is 

a  suffio 

lent  but  not  necessary 

condition  f < 

ii  the  ; 

solvability  <  :■ 

the  inverse 

problem . 
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Without  going  into  details  we  cite  the  example  given  in  Ref.  17.  For 
simplicity  taking  the  time  dependence  to  be  of  the  form  he(t)  =  &(t),  h^ft) 
=  6'(t),  where  '  =  d/dt,  and  assuming  the  field  to  be  of  the  form 

f+(x,t)  =  sin(k(x  -  ct) 3 ,  -a  <  x  -  ct  <  a,  ka  =  nn, 


the  source  function  is  found  to  be  given  by 

p(x,t)  =  -[k  &(t)  cos(kx)  -  c_^8'(t)  sin(kx)],  -a  <  x  <  a. 

Note  that  in  this  case  the  source  is  confined  to  a  finite  region  of  space, 

which  is  an  important  attribute  for  any  physically  realizable  source 

function.  {We  should  point  out  that  the  source  function  given  by  Moses  is 

erroneous  due  to  a  sign  error  in  evaluating  the  Fourier  transform  of  h  (t) 

o 

[his  Eq.  (2. 36')].} 

i i )  Solu t ion  of  the  Th ree-Dimensional  Wave  Equation  in  the  Wave  Zone 

The  general  solution  of  the  source-free  three-dimensional  wave 

.  .18 
equation  is 


f(r.t)-  ( 2 re) 


-3/2 


?  V 


ik. r  -iockt 


F(k , 


o  )  , 


(29) 


where  a  -  +! ,  and  k  =  jkj.  Making  use  of  the  method  of  stationary  phase, 
it  can  be  shown  that 


i  k .  r 

3  iru  e 
r 


_2  fti _  f  ikr 

kt  si’iG  L  e 


«r)M*  iK) 


ikr 


6(  0 


0’ 

X 


)  fi(  ♦ 


) 

A 


(30) 


Will! 

r  v  ( 0 

X 

’  V 

are  the  angles  in  the 

t  h  e 

a  tig  i.  t 

1 5>  1  n 

the  cone  along  /  ax i .*• 

cot  e  along  +z  axis,  and  O',  <t>  '  are 

x  x 

•j  i 

We  .shall  def  ine  cone*-  more 
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precisely  in  the  following.  Upon  substituting  Sq.  (30)  into  Eq .  (29),  the 
field  in  the  wave  zone  is  found  to  be  given  by 

,  1  rm  r  ik.(r-ct)_  ik(r+ct)  , 

f(r,  i)  =  (8n)  r  Im  j  dk  |e  f(k,6x,^x)  -  e  F(k.,0^  ,^)J. 

(31) 

Remarkably,  Eq.  (31)  shows  that  the  general  solution  of  the  three- 
dimensional  wave  equation  in  the  wave  zone  is,  apart  from  the  factor  1/r,  a 
superposition  of  one  dimensional  wave  motion,  expressed  as  functions  of 
f.  -  ct  and  r  +  et.  Moreover,  if  initially  f(r,t)  is  confined  to  a  given 
solid  angle,  ?(k,6,<f>)  will  be  significant  for  k  lying  within  that  angle 
and,  from  Eq.  (31),  the  solution  in  the  wave  zone  will  also  be  confined  to 
the  same  angle.  For  the  purposes  of  interpretation  it  is  convenient  to 
consider  propagation  in  cones,  as  indicated  in  Fig.  6.  In  particular,  an 
electromagnetic,  field  in  a  cone  is  defined  as  one  that  is  completely 
confined  to  a  cone  in  the  wave  zone. 

i i i )  Exact  So 1 u tion  of  t he_  Wave  Equation  from  the  Solution  in  the 
Wave  Zone  (Radon  Transforms ) 

The  purpose  of  this  subsection  is  to  point  out  that,  given  the 

solu tii  of  the  wave  equation  in  the  wave- zone  region,  it  is  possible  to 
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determine  the  solution  everywhere.  ’  In  particular,  one  seeks  the 

initial  conditions  f(r,t)  and  (S/St)f(r,t)  at  t  =  T.  This  is  referred  to 

as  the  inverse  initial  value  problem. 

The  solution  of  the  inverse  initial  value  problem  is  discussed  in  Kef. 

18  for  a  particularly  simple  case.  A  systematic  treatment  of  the  general 

i  9 

problem  is  possible  by  using  Radon  transforms. 

The  Radon  transform  F(l:,n)  of  a  function  t(r)  is  obtained  by 
integrating  f  over  all  planes  r.n  const. 
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F(k,n)  =  Jd^r  *(£)  S(r-n  -  k), 


A  24 

where  n  is  a  unit  vector.  The  usual  Radon  transform  is  defined  as  an 

integral  over  planes  whose  normals  vary  over  a  unit  sphere.  In  general, 

the  function  f(r)  defines  some  internal  distribution  (such  as  density)  of 

A 

an  object  and  F(k,n)  is  the  projected  distribution,  or  the  profile,  of  the 
object  on  the  plane  r.n  =  constant.  The  Radon  transform  is  a  very  useful 
tool  in  image  reconstruction  from  projections,  with  applications  in 
computer-assisted  tomography  (CAT)-scan,  radio  astronomy,  remote  sensing, 
etc . 

A  refinement  of  the  usual  definition  of  the  Radon  transform  shows  that 
only  the  transform  over  a  hemisphere ,  which  may  consist  of  disjointed 
parts,  is  sufficient  to  reconstruct  the  original  function.  It  can  then  be 
shown  that  the  task  of  obtaining  an  exact  solution  of  the  three-dimensional 
wave  equation  from  the  solution  in  the  wave  zone  reduces  to  tak'ng  the 
inverse  of  the  refined  Radon  transform  of  the  solution  in  the  wave  zone. 

From  Eq.  (31)  it  is  known  that  in  the  wave  zone  the  solution  of  the 
wave  equation  is,  apart  from  a  factor  1/r,  a  function  of  only  r-ct  or  r+ct, 
representing  propagation  along  rays  confined  to  a  cone.  The  field  in  the 
wave  zone,  therefore,  defines  the  range  of  the  unit  vector  n  and  the 
amplitude.  This  information  is  essentially  equivalent  to  knowing  the 
projections  in  different  directions.  The  exact  field  may  then  be 
reconstructed  from  the  set  of  projections, 
iv)  Example  of  a  Bullet 

We  close  Section  VI 1  by  discussing  an  explicit  example  of  a  bullet 

2  ( ) 

which  is  a  solution  of  the  scalar  wave  equation, 
electromagnetic  bullet  is  given  in  Ref.  20. 


An  example  of  an 


A  bullet  which  is  confined  to  a  finite  volume  in  the  wave-zone  is 


given  by 

f(r,t)  =  -  0)  [  h(  r-a-ct )  -  h(  r-b-ct )  ] /r ,  r  -*  ®,  (32) 


where  h  is  the  Heaviside  function,  a  and  b,  with  b  >  a,  denote  the 
boundaries  of  the  bullet  along  the  radius,  and  2 a  is  the  vertex  angle  of 
the  cone  containing  the  bullet-  Note  that  this  solution  is  causal  and  of 
finite  energy,  and  has  a  form  which  is  consistent  with  the  remarks 
following  Eq.  (31).  It  represents  a  packet  of  energy  "shot"  through  a  cone 
whose  axis  coincides  with  the  z  axis. 

This  solution  can  be  easily  verified  by  letting  f  =  g(r,9,t)/r  and 
noting  that 
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Assuming  g(r,0,t)  =  h(r-ct-ro)Y(0) ,  one  finds 
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a^_ 

at' 


r^sin© 
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(sin©  ^  Y 
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(33) 


which  tends  to  zero  for  all  continuous  Y(0)  as  r  -+  For  the  case  when 

2  2  1/2 

Y  ( 0)  =  n(  a  -  0),  we  let  dY/d©  =  exp[-(a-©)  /2A  )/(2fl)  A,  with  A  -»  0. 

Note  that  as  A  ->  0,  dY/d©  -»  -S(o-0)  =  dn(a-0)/d0.  Upon  substituting  this 
form  into  Eq.  (33),  we  notice  that  as  A  is  made  to  approach  zero,  r  must 
increase  indefinitely  in  order  for  the  right-hand  side  of  Eq .  (33)  to 
approach  zero.  This  shows  very  clearly  that  Eq .  (3?)  is  indeed  a  solution 
in  the  wave  zone. 

The  spot  size  may  he  defined  by: 

o  ■>  n 

i‘  Jdid0d4>  f  i ‘sin©  (2rsin©)‘  /  J  d  r  d  0d  4>  f  i‘sin0,  (34) 
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where  2rsin0  is  the  width  of  a  cone  of  half -angle  0.  Substituting  Eq.  (32) 
into  Eq.  (34),  one  obtains 

r s  -a>  [  ( 8  -  9co so  +  cos3cr)/3(l  -  eoscr)  |  ^^c t ,  t  ->  <»,.  (35) 

indicating  a  linear  increase  with  time  for  the  spot  size,  as  in  the  case  of 
the  Gaussian  pulse  in  Eq.  (14b),  although  the  constants  of  proportionality 
are  different.  We  mention  in  passing  that  the  solutions  given  by  Moses  and 
Prosser  are  distinguished  from  the  other  solutions  reviewed  here  by  not 
having  an  explicit  dependence  on  the  frequency  or  the  wave  number. 

As  mentioned  in  the  previous  subsection,  to  obtain  the  exact  solution 
everywhere  one  has  merely  to  evaluate  the  inverse  Radon  transform  of  the 
solution  in  the  wave  zone,  Eq.  (32).  Since  the  derivation  is  lengthy  we 
shall  simply  quote  the  result  The  exact  solution  is  given  by 

f(r,t)  f a( r ,  t )  -  fb(r,t),  (36a) 

where 

f  (r,t)  =  h(a-0)  {h[a+ct-rcos(0-cr)  ]  -  f|(a  +  c  t-r ) } /r 
3  — 

+  v  {b[a+ct-rcos(6+cr)  j  -  fj[ a+ct-rcos( 0- o)  ] )  / nr ,  (36b) 

3 


wi  th 


v  =  cos  ^((cosa  -  cosg  cos0)/sing  sin0],  0  <  v  <  n, 

3  3  3  3 

and 

g  -  cos  J|(a+ct)/r|,  0  <  g  <  n/2, 

3  cl 

and  where  vj(r, t)  is  identical  to  t  (r,t)  except  that  a  is  replaced  by  b. 
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The  field  in  Eq .  (36)  is  identically  zero  for  all  r  downstream  of  the 
bullet,  i.e.,  for  all  r  <  a  +  ct.  The  wave-zone  limit  is  obtained  by 
taking  r,  t.  ->  ”.  Then,  since  Eq .  (36)  corresponds  to  propagation  in  the 
positive  cone  (Fig.  6),  taking  r  -  ct  =  constant  Eq .  (32)  is  recovered. 

The  requirement  r  -  ct  =  constant  is  equivalent  to  observing  the  bullet  in 
a  co-moving  frame.  Close  to  the  origin  the  exact  solution  in  Eq .  (36) 
spreads  out  of  the  cone  significantly.  However,  in  the  wave  zone  the 
solution  is  confined  to  the  cone  and  is  independent  of  the  angle  0  therein. 
Finally  it  has  been  shown  in  Ref.  20  that  the  difference  between  the  exact 
solution  in  Eq.  (36)  and  the  wave-zone  solution  in  Eq .  (32)  becomes  small 
quite  rapidly  as  r,  ct  increase,  with  r  -  ct  held  fixed. 

In  principle,  one  can  now  use  the  inverse  source  method  to  determine 
the  sources  that  lead  to  the  bullet  described  by  Eq.  (32).  To  our 
knowledge,  however,  this  computation  has  yet  to  be  performed. 


VIII.  Summary  and  Concluding  Remarks 

The  motivation  for  much  of  the  research  reviewed  herein  stems  from 
the  need  to  propagate  a  beam  of  radiation  over  long  distances  without  an 
appreciable  decrease  in  the  intensity.  Possible  applications  would 
include:  power  beaming,  advanced  radar,  laser  acceleration  of  particles 

and  directed  energy  sources.  This  need  has  led  to  a  great  deal  of  interest 
in  such  fundamental  subjects  as  diffraction  and  new  solutions  of  the  wave 
equation. 

It  has  been  reiterated  that  the  physical  basis  for  diffraction  of 
waves  is  the  well-known  relation  Ak^Ax^  =•-  1,  for  i  =  1,2,3.  By  virtue  of 
this,  it  is  simple  to  determine  the  scale  length  for  the  diffractive 
spreading  of  a  beam  with  an  arbitrary  transverse  profile.  Thus,  a 
knowledge  of  the  spectrum  is  sufficient  to  determine  the  maximum 
propagation  distance  of  the  beam.  Since  diffraction  is  unavoidable,  by 
concentrating  the  energy  in  the  high  frequencies  one  can  only  delay  the 
spreading  of  the  beam. 

Four  examples  of  the  research  effort  on  the  subject  of  beam 
propagation  have  been  reviewed  herein.  The  conclusions  are  as  follows. 

i)  Electromagnetic  Missiles 

Experiment  indicates  that  a  suitably  tailored  pulse  shape  can  be 
designed  to  have  an  energy  decay  rate  limited  by  the  highest  frequencies 
present  in  the  pulse.  This  is  fully  consistent  with  our  understanding  of 
d i f  f rac t  i on . 

ii)  Bessel  Beams 

It  is  .shown  that  as  fat  as  propagation  is  roiieei  ned  Bosse!  beams  are 
not  "resistant  to  the  diffractive  spreading  commonly  associated  with  all 
wave  ptopagat ion".  These  beams  propagate  no  luithet  that  Gaussian  beams  or 
plane  waves  with  the  same  transverse  d imeus ions . 


i i i )  Elec tromagnet  ic  Pit: ected  Energy  Pulse  Trains 

The  diffractive  properties  of  the  pulse  form  studied  most  intensively 
under  this  general  heading  are  described  by  diffraction  theory.  These 
pulse  trains  do  not  "defeat  diffraction". 

i v )  Elec tromagnetic  Bulla ts 

Given  a  radiation  wave  packet  in  the  wave-xone  which  is  confined  to  a 
suitable  solid  angle  and  extends  over  a  finite  radial  extent,  one  can 
determine  the  source  required  to  generate  the  wave  packet.  As  of  this 
writing,  however,  this  problem  has  not  been  solved  for  a  practical  case. 


Acknowledgment 

This  work  was  supported  by  the  Office  cf  Naval  Research  and  the 
Department  of  Energy. 


Appendix 


The  purpose  of  this  appendix  is  to  examine  the  transition  from  the 

Fresnel  to  the  Fraunhofer  region  for  a  clipped  Bessel  beam  and  a  clipped 
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Gaussian  beam  within  the  context  of  the  Huygens-Fresnel  approximation. 

The  clipping  is  assumed  to  be  caused  by  a  finite-size  aperture. 

In  the  case  of  the  Bessel  beam  the  field  distribution  at  the  aperture 
has  the  form  u(  r , z-0)  =  Jo(k^r)  within  a  circular  aperture  of  radius  d. 
Making  use  of  Eq.  (6)  the  amplitude  at  a  point  on  the  axis  of  symmetry  is 
given  by 

d  2 

«  2  dr'  r'  JQ{k^r' )  22,6 .  (Al) 

o 

The  limit  k^  =  0  corresponds  to  the  case  of  plane  waves  incident  on  a 
circular  aperture,  as  in  Section  III.  The  intensity,  given  by  Eq .  (9), 
fails  off  rnonotonically  for  z  >  Zp,  where 

Zp  =  2  rid  2  /  X .  (A2 ) 

For  kj  >  0,  the  oscillatory  behavior  of  the  Bessel  function  in  Eq .  (Al) 
tends  to  phase  mix  the  integrand,  effectively  reducing  the  upper  limit  of 
the  integration.  Consequently,  the  boundary  of  the  Fresnel  region,  beyond 
which  the  radiation  appears  to  be  emitted  essentially  from  a  point  source, 
is  reached  prior  to  Zp  =  2nd"/A. 

Since  there  is  no  simple  analytical  appt  ox i mat  ion  to  the  integral  in 
Fq .  (Al),  consider  the  case  of  a  cosine  beam  u(x,y,z  0)  cos ( k  x  )  cos ( k  y ) , 
which  i s  the  cartesian  equivalent  of  a  Bessel  beam.  The  a pet  t  ut  e  is  a 
i octangular  opening  in  the  x  y  |>latte  detined  hy  \  >'■  |  <  X,  |y|  <  Y). 


IS 


In  terms  of 


(A3) 


l+  =  (k/2z)1/2(X  ±  kxz/k>, 


the  amplitude  on  the  2  axis  has  the  form  \fi  =  I  I  ,  where 

x  y’ 


_ 2  2  /  2  k. 

I  =  e  X  [C(l  )  +  iS(E  )  +  ca  )  +  isa  )], 

A  +■+■—  — 


(A4) 


and 


t 

1/2  f  ? 

S(t)  =  ( 2 / it)  dt  sint  , 

o 


t 

C(t)  =  (2/rt)1/2J  dt  cost2, 
o 


are  the  Fresnel  integrals.  The  expression  for  I  is  obtained  from  that 
for  1^  by  making  the  replacements  ->  k^,  X  -*  Y. 

We  are  interested  in  the  intensity  well  within  a  nominal  Fresnel 
region  defined  by  the  width  of  the  aperture 


z  <<  k(X2,Y2)/2, 


(A5) 


but  sufficiently  far  from  the  aperture  so  that  the  radiation  diffracted 
from  one  edge  can  reach  the  z  axis: 

z  -  k(X/kx,Y/k  >.  (A6) 


Taking  the  appropriate  limits  of  the  Fresnel  integrals,  the  intensity  is 
I  =  I  /16,  where  1  is  the  intensity  at  the  nit  liar  ting  aperture.  This 
analysis  indicates  that  Eq.  (Ab)  defines  the  boundary  ol  the  true  Fresnel 
region.  The  terms  proportional  to  1  k  z/k  in  F.q .  (A3)  represent 
propagation  at  an  angle  ■>  sin  *  (k  /k)  to  the  ::  axis,.  As  a  consequence  the 
drop  in  intensity  characterising  the  transition  to  the  Kiaimhoiti  region 
takes  plane  at  the  location  indicated  by  F.q.  (AO)  rattier  than  by  tin-  1  iglrt 
hand  side  of  F.q.  (Ah). 


If. 


Returning  to  the  Bessel  beam  with  d  >>  1/k.j ,  in  the  region  of 
significant  phase  mixing  in  the  integrand  of  Eq.  (Al)  the  Bessel  function 
has  the  asymptotic,  form  J  (z)  ~  (2/nz)  "cos(z  -  tt/4).  Substituting  this 
into  Eq.  (Al)  and  comparing  the  phases,  it  follows  that  the  Fraunhofer 
region  for  the  Bessel  beam  commences  at  ZR  =  2 nd / X .  From  the  definition 
of  Zp  in  Eq.  (A2)  we  note  that  ZR/Zp  =  l/'k^d  <<  1  for  the  experimental 
parameters  in  Ref.  11.  Thus,  we  see  that  a  Bessel  beam  is  not  optimum  as 
far  as  the  diffractive  fall-off  of  the  intensity  is  concerned. 

For  the  Gaussian  beam,  substituting  u( r)exp(-r2/w2 )  into  Eq.  (6)  and 
performing  the  integral,  the  intensity  on  the  axis  of  symmetry  is  found  tc 
be  given  by 

-2(d/w  )2  -<d/w  )2  - 

x  a  1  r  e _ -  2e _ cos(k.d  /2z) ,  (A7) 

1  +  (2/Zr)2 

where  Z  is  the  Rayleigh  range  defined  in  Eq.  (1).  For  the  parameters  in 

Ref.  11,  d/wo  >>  1,  and  the  scale- length  for  the  intensity  to  drop  to  a 

quarter  of  its  initial  value  is  on  the  order  of  2ZR,  as  in  Eq .  (4c).  This 

same  scale  length  is  roughly  applicable  to  the  case  of  a  wider  Gaussian 

beam  with  w  =  d.  For  an  infinitely  wide  beam,  w  -+  and  Eq .  (A7)  goes 
o  o 

over  to  the  case  of  plane  waves,  Eq.  (c>). 

Fot  the  off  axis  intensity  of  the  cosine  beam  we  limit  the  discussion 
to  the  case  where  the  observation  point  is  an  integral  number  of  half 
periods  oil  of  the  z  axis.  In  analogy  with  Eq .  (AT)  we  define 

K  (  k / ?  <’.)'  7 '  |  X  >  (1:  :  ■  k  n  n/k  1  i  , 

+  >:  >- 

n  a.  ?.  )!  ’  i  x  .  av.-  i  •  i:  n  i  v  >  i . 


where  n  ,Jl/k  ,  with  n  an  integer,  is  the  x-coordinate  of  the  observation 

XX  X 

point.  We  assume  n  Jl/k  <  X.  In  terms  of  these  variables,  the  amplitude 

has  the  form  K  I  I  ,  where 

x  y 


I  «  |C(F  )  +  iS(E  )  +  C(f  )  +  iS(f.  ) 


+  c<n  )  +  isov)  +  c(n  )  +  is(n  )], 


(AS) 


where  C  and  S  are  the  Fresnel  integrals  defined  earlier  and  the  expression 

for  I  is  obtained  from  tnat  for  I  by  making  the  substitution  k  -»  k  ,  X  -» 
y  x  J  &  x  y  ’ 

Y,  and  n  -*  n  . 

x  y 

For  plane  waves  (k  -»  0)  and  on  the  symmetry  axis  (n  ->  0)  the 

x , y  ’  '  x ,  y  ' 

transition  from  the  oscillatory  to  the  monotonically  falling  behavior  of 

the  Fresne1  integrals  in  Eq .  (A8)  taker  place  at  z  «  k(X^,Y^)/2.  This 

marks  the  boundary  between  the  Fresnel  and  Fraunhofe*  regions.  For  the 

cosine  beam  and  on  the  symmetry  axis  (n _  ->  0)  Eq .  (A8)  reduces  to  Eq. 

x ,  y 

(A4)  and  hence  Eq .  (A6)  defines  the  boundary  between  the  tvr  regions.  Away 
from  the  symmetry  axis  (n  >  0)  the  behavior  is  somewhat  more 

x  •  y 

complicated.  On  the  aperture,  z  =  0,  and  F  ,  »i+  ->  ®,  whence  I  =  2(1  +  i). 
For  z  sufficiently  large  so  that  0  <  t|  <1,  but  F,  ,  F  ,  0  >>  1,  the  last 

—  -f  —  + 

two  terms  in  Eq.  (A8)  are  small  and  1  -»  3(1  +  i)/2.  From  the  definition 
of  n  ve  note  that  as  the  bservation  point  approaches  the  edge  of  the 


aperture  (n  n/k  -»  X)  this  reduction  in  the  value  of  j  is  obtained  at 

XX  X 

smaller  values  of  z,  according  to  the  formula  z  --  k(X  -  n  tl/k  )k  .  Figure 

,2 


/,  Wf; ;  r  h  plot 
parameters  similar  to  that  iri  Re 


1  I'  as  a  f  unc  tian  of  z,  con f  '  t  nus 
x  1 


1  1 


this  behavior  for 
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